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ABSTRACT 

The velocity divergence power spectrum is a key ingredient in modelling redshift space 
distortion effects on quasi-linear and nonlinear scales. We present an improved model 
for the z — velocity divergence auto and cross power spectrum which was originally 
suggested by Jennings ct al. 2011. Using numerical simulations we measure the velocity 
fields using a Delaunay tessclation and obtain an accurate prediction of the velocity 
divergence power spectrum on scales k < 1/iMpc -1 . We use this to update the model 
which is now accurate to 2% for both Pee and Pes at z = on scales k < 0.65/iMpc -1 
and k < O.^bhMpc" 1 respectively We find that the formula for the redshift dependence 
of the velocity divergence power spectra proposed by Jennings et al. 2011 recovers 
the measured z > P(k) to markedly greater accuracy with the new model. The 
nonlinear Pgg and Pes at z = 1 are recovered accurately to better than 2% on scales 
k < 0.2/iMpc _1 . Recently it was shown that the velocity field shows larger differences 
between modified gravity cosmologies and ACDM compared to the matter field. An 
accurate model for the velocity divergence power spectrum, such as the one presented 
here, is a valuable tool for analysing redshift space distortion effects in future galaxy 
surveys and for constraining deviations from general relativity. 

Key words: Methods: iV-body simulations - Cosmology: theory - large-scale struc- 
ture of the Universe 



1 INTRODUCTION 

In the hierarchical model of structure formation, gravita- 
tional collapse and the accelerating cosmic expansion are 
two competing effects which determine the rate at which 
structures, such as galaxies and clusters, grow in the Uni- 
verse. In addition to the Hubble flow, galaxies possess pe- 
culiar velocities, arising from inhomogeneities in the lo- 
cal density field, which can be used t o probe the growth 



rate of structure JPercival et aZ.1 12007 : IGuzzo et al. 120081: 



Blake et al] |2010| . 1201 ll : ISchlegel et al\W0l\ ; iBeutler et al\ 
20121 ). These peculiar velocities distort the clustering sig- 
nal along the line of sight giving rise to redsh ift space 



distortions (see e.g. iKaiserl 1 19871 : iHamiltonl 1 19981 ). It has 
been shown that a key ingredient to improve models of 
the power spectrum in redshift space is the inc lusion of the 
nonlinear velocity divergence power spe ctrum (|Scoccimarrol 
12004 IJennings. Baugh fc Pascolill2011al lbl) . In this paper we 
present an improvement to the model for both the auto 
and cross velocity divergence power spectrum presented in 



E-mail : cj ennings® kicp . uchicago . cdu 



IJennings. Baugh fc Pascolil (|2011al ). The refinement of the 
model is driven by using a volume weighted Delaunay tesse- 
lation method to measure the nonlinear velocity fields to 
greater accuracy and to smaller scales then presented in 
IJennings. Baugh fc Pascolil (|2011al ). 

Measurements of the growth rate of structure can be 
used to determine if the accelerating expansion is the re- 
sult of a dark energy component which behaves as a repul- 
sive form of gravity or if Einstein's t heory of gravity breaks 
down on cosmological scales (see e.g. iBertschinger fc Zukinl 
120081 ). Independent measurements of the growth rate can 
be obtained by measuring the clustering of galaxies in 
redshift space and recently there has been renewed in- 
terest in improving the models for the clustering sig- 
nal in redshift space (|Scoccimarrd 2004: iPercival fc White! 
20091: iTaruva et al.ll2010l: I Jennings. Baugh fc Pascolill2011al: 
Seliak fc McDonald! l201ll; iTang et o.Z.1 l201ll: iReid fc White! 



201ll : lKwan. Lewis fc Lmderl2012l )7 Also recently it has been 



pointed out that the velocity divergence power spectrum is a 
more sensitive probe of modified gravity than the nonlinear 
matter power spectrum (| Jennings et all |2012| ) but models 
for the 2D redshift space power spectrum are currently not 
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accurate enough to allow this quantity to be extracted to 
a high precision. However, given an accurate model for the 
velocity divergence power spectrum in a ACDM cosmology, 
this can be used to construct predictions for the moments of 
the power spectrum in redshift space which can be directly 
compared with data and has been shown to be accurate 
in detect ing deviations from general relativity using sim- 
ulatio ns l|jennings. Baugh fc Pascoll l2011bl ; Ijennings et all 
12012ft . 

The linear continuity equation, V • v = —aHfS, gives a 
one to one correspondence between the velocity and density 
fields where the linear growth rate / = dln5/dlna is a cos- 
mology dependent factor, 8 is the matter overdensity, v is the 
peculiar velocity and H is the Hubble parameter. Once over- 
densities become no nlinear, this relationship no longer holds. 
iBernardeaul (|l992T ) derived the nonlinear relation between 8 
and V ■ v in the case of an initially Gaussian field. Many au- 
thors since then have extended this relation using e.g. higher 



2.1 Simulation details 



e.g. IScoccimarro & Friemanl 


19961; 1 Chodorowski & Lokas! 


199tT Kudlicki et al.l 20001; 


Bilicki & Chodorowskil 20081; 


Kitaura et al.l l201lh. ICiecielg 


& Chodorowski 1 ( 2004T) were 



the first to apply numerical simulations to model the non- 
linear velocity and density spectra, instead of using pertur- 
bation theory. Using a grid-based pressureless hydrodynamic 
code their simulations had limited resolution but were able 
to fit the ratios up to k = 1/iMpc -1 and show that pertur- 
bative predictions fail at k > 0.2/iMpc . 

The density veloc ity relation considered in this pape r 
was first presented in Ijennings. Baugh fc Pascoll (|2011al ). 
where the parameters were obtained by fitting to a mass 
weighted measurement of the velocity divergence field. In 
this paper we present updated parameters obtained by 
fitting to the volume weighted velocity divergence power 
spectra which have been measured accurately on smaller 
scales than was possible with the mass weighted method 
(IPueblas fc Scoccimarr oll2009l). We also show that the model 
proposed bv lJennings. Baugh fc Pascolil (|2011al ) for the red- 
shift dependence of the density velocity relation recovers the 
velocity power spectra at z > with greater precision using 
the new parameters. 

This paper is organised as follows: In Section 12.11 we 
describe the iV-body simulations used in this paper. In Sec- 
tion [221 we discuss two methods for measuring the velocity 
divergence power spectrum. In Section|3]we present the main 
results of this paper including the updated parameters for 
the nonlinear velocity divergence power spectra. 



2 MEASURING P ge AND P gs FROM N-BODY 
SIMULATIONS 

In section 12.11 we present the details of the N-body simula- 
tions carried out and discuss several methods which can be 
used to measure the velocity divergence power spectrum in 
Section 12.21 We focus on a Delaunay tesselation approach 
and show the accuracy which with this method can recover 
Pge and Pes- 



We u se the N-body simulations carried out by iLi et all 
These simulations were performed at the Institute of 
Computational Cosmology using a modified vers ion of the 
mesh-based iV-body code RAMSES (|Tevssierl [20021 1. Assum- 
ing a ACDM cosmology, the following cosmological param- 
eters were used in the simulations: fi m = 0.24, SIde = 0.76, 
h = 0.73 and a spectral t ilt of ns = 0.961 (in agreement with 
e.g. ISanchez et fflZ~ll2009r i. The linear theory rms fluctuation 
in spheres of radius 8 /i -1 Mpc is set to be as = 0.769. 

The simulations use N = 1024 3 dark matter particles to 
represent the matter distribution in a computational box of 
comoving length 1500fe~ 1 Mp c. The ini tial conditions were 
generated at z = 49 using the MPgrafic code. Note the non- 
linear matter power spectrum is measured from the simula- 
tions by assigning the particles to a mesh using the cloud in 
cell (CIC) assignment scheme and performing a fast Fourier 
transform (FFT) of the density field. To compensate for the 
mass assignment sch eme we perform an ap proximate de- 
convolution following iBaumgart fc Frvl j 199 if ). 

Throughout this paper the velocity divergence is nor- 
malized to 9 — — V • v/(aH /), where v is the peculiar veloc- 
ity, / = dln<5/dlna is the linear growth rate, H is the Hubble 
parameter and a is the scale factor. Using this normaliza- 
tion 8 is dimensionless and the linear continuity equation is 
given by 6 = 8. 



2.2 Measuring the velocity divergence field 



Measuring the 
from numerical 
difficult if 



Scoccimarrcl (120041') 



velocity divergence field accurately 
simulations on small scales can be 
we ighted approach is used 



Jennings. Baugh 



.Pueblas fc Scoccimarrcl (|2009l ); 
Pascolil (|2011al . for example) . Some 



volume weighted measures of the velocity field have also 
been de veloped (see e.g. Bernardeau fc van de Weygaerj 
1 19961 ; IColombi. Chodorowski fc Tevssierj l2007f ) in- 
cluding the Delaunay te ssellation field estima- 
tor (DTFE) method dSchaap fc van de WevgaertJ 
20001; Ivan de Weygaert fc Schaad 120091 ; 



Cantun fc van de Weygaerj 2011 ) 



In the mass weighted approach, si mply interpola t ing th e 
velocities to a grid, as suggested by IScoccimarr 3 1120041) . 
gives the momentum field which is then Fourier trans- 
formed and divided by the Fourier transform of the den- 
sity field, which results in a mass weighted velocity field 
on the grid. One of the main problems with this approach 
is that the velocity field is artifically set to zero in regions 
where there are no particles, as the density is zero in these 
empty cells. IPueblas fc Scoccimarrcl (|2009T l also found that 
this method does not accurately recover the input velocity 
divergence power spectrum on scales k > 0.2/iMpc" 1 in- 
terpolating the velocities of 640 3 particles to a 200 3 grid. 
Using simulations of 1024 3 particl es in a 1.5/i -1 Gpc box, 
Ijennings. Baugh fc Pascolil (|2011af ) found that the maxi- 
mum grid size that could be used was 350 3 without reaching 
the limit of empty cells. 

The limit on the maximum grid size which can be used 
in a mass weighted estimate o f 6 means that the veloc- 
ity dive rgence P{k) presented in ljennings. Baugh fc Pascolil 
(|2011al ) is only accurate on large scales k < 0.2/iMpc _1 . 
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Figure 1. Left panel: Ratios of the 2 = nonlinear matter, Pgg, and velocity divergence, Pgg, power spectra to the matter power 
spectrum at 2 = 4, P$$(z = 4), scaled using the ratio of the square of the linear growth factor at z = 4 and z = for ACDM, measured 
from six simulations are shown as filled blue circles and empty green diamonds respectively. The ratio of the Pgg/Pss measured from the 
simulations at z = is shown as empty purple squares. The shaded grey region shows the errors on this ratio measured from the scatter 
amongst six simulations. Right panel: Similar to the left panel. The ratios Pgs{z = 0)/Pss( z = 4) and Pgs(z = 0)/Pss(z = 0) are shown 
as empty green diamonds and purple squares respectively. Note 8 = — V ■ v/(aHf). 




Figure 2. Left panel: The ratio of the fitting formula given in Eq. [T]for the velocity divergence P(k) to the measured power spectrum, 
Pgg/Pgg m , at z = (green solid), z = 0.4 (red dashed) and 2 = 1 (blue dot dashed). Right panel: Similar results as shown in the left 
panel but for the cross power spectrum Pgg. The horizontal dotted grey lines show a region of 2% accuracy in the fit. 



Also, it is not clear on what scales the discrepancy be- 
tween the actual and measured ve locity divergence power 
spectr um, which was found by IPueblas fc Scoccimarrol 
l|2009h on scales k < 0.02/iMpc" 1 (see Figure 13 
that p aper), is present in the ve locity P(k) mea- 



sured by Jennings. Baugh fc Pascolil (|20Tlal ). However, 



Ijennings. Baugh fc Pascolil (|2011al ) find that the measured 
Pgg and Pg$ agree with the linear continuity equation on 
large scales within the errors from eight simulations (see 
Figure 3 in that work). 

The limitations of using the mass weighted method to 
measure Pgg implies that t he fitt ing formula presented in 
Ijennings. Baugh fc Pascolil (|2011al ) is only valid over a lim- 
ited range of scales. The main aim of this paper is to present 
an improved version of the formula for Pgg and Pgs given by 
Ijennings. Baugh fc Pascolil (|2011al ) by fitting to power spec- 
tra which have been measured using the DTFE method This 
code constructs the Delaunay tessellation from a discrete set 
of points and interpolates the field values onto a user defined 



grid. For the L box = 1500/i -1 Mpc simulation we generate 
the velocity auto, Pgg, and cross power spectrum, P$g, on 
a 1024 3 grid. The velocity divergence field is interpolated 
onto the grid by randomly sampling the field values at a 
given number of sample points within the Delaunay cells 
and than taking the average of those values. The resolution 
of the mesh used in this study means that mass assignment 
effects are negligible on the scales of interest here. 

In Fig. [l]we show the ratios of Pss(z = 0)/Pss(z = 4) 
and Pgg(z = 0)/Pss(z = 4) which represent the average 
ratio measured from six simulations as filled blue circles 
and empty green diamonds respectively. Here Pss(z = 4) 
has been scaled using the ratio of the square of the lin- 
ear growth factor at z = 4 and 2 = for ACDM. 
Plotti ng the ratio in this way r emoves the sampling vari- 
ance (|Baugh fc Efstathioul Il994l ) and shows the agreement 
of the measured velocity and matter P(k) with the pre- 
dictions of the linear continuity equation on large scales, 
k < 0.02/iMpc- 1 for Pg 9 (z = 0)/P ss (z = 4) and k < 
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0.09/iMpc _1 for P ss (z = 0)/P S s(z = 4). This result agrees 
with similarly measured rat i os fro m simulations for the mat- 
ter (see e.g. lAngulo et al\ l2008h and the velocity diver- 
gence power spectrum Ijennings, Baugh fc Pascoll l2011al ; 
IJennings er^Z|]2012l : [Li et al. Il2012t ). 

The ratio of Pee(z = 0)/Pss(z = 0) is plotted in the 
left panel of Fig. [T] as empty purple circles. The grey shaded 
region shows the errors measured from the scatter amongst 
six simulations. In the right panel of Fig. [1] we plot similar 
ratios as in the left panel but for the velocity divergence 
cross power spectrum Pes. Note the increase in the ratio 
of Pes(z = 0)/Pss(z = 4) on scales fc > 0.1/iMpc" 1 which 
is due to the nonlinear growth in the matter field on these 
scales. Overall the results plotted in Fig. [T] show that the 
velocity divergence field computed using the DTFE method 
agrees with the predictions of linear theory on extremely 
large scales and allows us to accurately measure the power 
spectra on nonlinear scales provided we use a large enough 
grid to interpolate the velocities. We have checked that the 
effect on the measured velocity power spectra using smaller 
grid sizes such as 256 3 and 512 3 and find excellent agreement 
between the P(k) measured using a 512 3 and 1024 3 grid on 
scales k < 1/iMpc -1 ; the P(k) measured using a 256 3 grid 
shows deviations at fc ~ 0.7/lMpc _ . 



3 RESULTS 

Using the measured nonlinear velocity divergence auto and 
cross power spectrum on a 1024 3 grid presented in Section 
I2.2l and the nonlinear matter power spectrum measured from 
the simulations at z = 0, we can fit for new paramete rs in 
the following model (|jennings. Baugh fc Pascolill2011ah : 



P*v(k) = 9(Pss(k)) = 



aoVPss(k) + aiPss(k) 
«2 + a 3 P S s(k) 



(1) 



where Pss is the nonlinear matter power spectrum. We ob- 
tain the following: 

• Pxy = Pse- 

a = -12483.8, ai = 2.554, a 2 = 1381.29, q 3 = 2.540; 

• Pxy = Pee ■ 

a = -12480.5, Qi = 1.824, a 2 = 2165.87, a 3 = 1.796; 

all points were weighted equally in the fit. Note the pa- 
rameters ao-3 are not dimensio nless and in this work their 
units d iffer from those quoted in ljennings. Baugh fc Pascolil 
(|201 lal ~) as the best fit parameters presented here generate 
velocity divergence power spectra which are normalized as 
6 = V • v/(aHf). Note the growth rate at z — for the 
ACDM cosmology considered in this work is / = 0.452. We 
do not give the units of «o-3 here but these can be easily 
found given this normalisation for theta. The power spectra 
used for this fit are the average Pee, Pse and Pss measured 
from six ACDM simulations. 

In Fig. H we plot the ratio of Pfe/Pet and Pfs/Pef" 
in the left and right panels respectively at z — as a green 
solid line. Here the prefix 'fit' denotes the velocity diver- 
gence power spectrum found using the new best fit parame- 
ters and the nonlinear matter P(k) measured from the sim- 
ulations in Eq. [T] The prefix 'sim' denotes the nonlinear ve- 
locity divergence measured directly from the simulations us- 
ing the DTFE. The formula given in Eq.[T]is accurate to 2% 



(shown as the region enclosed by the dotted gray lines in Fig. 
[IJ on scales fc < 0.65/iMpc" 1 for Pee and fc < 0.35/iMpc _1 
for Pes at z = 0. Note the lower limit for the domain of 
the function given in Eq.EQis fc ^ 0.006/iMpc" 1 . Although 
the values for ao-3 were obtained by fitting to a ACDM 
simu lation with a particular se t of co smological parame- 
ters, IJennings. Baugh fc Pascolil |2011af ) found that this re- 
lation between the density and velocity power spectra is 
quite insensitive to both the cosmological model, for smooth 
dark energy models such as quintessence, and the choice 
of cosmological par ameters. The s e resu lts ag ree with previ- 
ous st udies such as iBernardeau I (|l992h and iBouchet et all 
(|l995h . 

The redshift dependence of Pee and Pes 
can be described using th e following formula 
(|jennings. Baugh fc Pascolill2011al ): 



Pxy{k, Z ) — 



g(P S s(k,z = 0))-P ss (k,z = 0) 
c 2 (z = 0,z') 

+Pss{k, z) , 



(2) 

where g(Pss) is the function in Eq. [T]and P xy is either the 
nonlinear cross or auto power spectrum, Pse or Pss and the 
function c is given by 

D(z) + D 2 (z) + D' A (z) 



c(z, z ) = 



(3) 



D(z')+D 2 (z') + D 3 {z') ' 

and D(z) is the linear growth factor. 

In Fig. [2] we show the ratio of Pee from Eqns. [T]&[2]to 
the measured velocity divergence power spectrum P|™ at 
z = 0.4 (red dashed) and z = 1 (blue dot dashed). All of 
the results shown in this plot represent the average over six 
ACDM simulations. It is clear from Fig. [2] that the model 
given in Eq. [5] is accurate to better than 2% on scales fc < 
0.2/iMpc" 1 and to 20% for 0.2 < fc^Mpc" 1 ) < 0.4 at z = 1. 



4 SUMMARY 

Measuring the growth of structure using the anisotropic 
clustering signal in redshift space is an important tool for 
discriminating between dynamical dark energy or modified 
gravity and the standard ACDM cosmological model. The 
velocity divergence power spectrum has been shown to be an 
important ingredient in modelling redshift space distor tions 
(|Scoccimarrd 120041 ; IJennings. Baugh fc Pascolil l2011al ) and 
shows larger deviations between modified gravity cosmolo- 
gies and ACDM than the differences found in the nonlinear 
matter power spectrum using numerical simulations (see e.g . 
Sch midt et al. I l2009bl ; IJennings ef~a/1l2012l ; iLi et al. II20I2I ) . 

In this paper we present an improved model for the z — 
velocity divergence auto and cross power spectrum which 
was o riginally suggested by IJennings. Baugh fc Pascolil 
|2011al ). By measuring the velocity fields using a Delaunay 
tesselation we obtain an accurate prediction of the velocity 
P(k) to fc ~ 1/iMpc -1 which we use to update the parame- 
ters for the fitting formula given in Eq. [1] We find that this 
model is accurate to 2% on scales fc < 0.65/iMpc _1 for Pee 
and fc < 0.35/iMpc" 1 for Pes at z = 0. We find that the 
formula for the redshift dependence of the velocity power 
spectra given in IJennings. Baugh fc Pascolil (|2011al ) recov- 
ers the z > P(k) to markedly greater accuracy with the 
new parameters. This model for the redshift evolution of Pee 
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and Pes, given in Eq.[2l is accurate to less than 2% on scales 
k < 0.2/iMpc -1 and to 20% for 0.2 < fc(/iMpc _1 ) < 0.4 at 

2=1. 

The improved model presented here accurately de- 
scribes the nonlinear density velocity relation and allows 
the velocity divergence power spectra to be easily and accu- 
rately predicted over a range of scales k < 0.9/iMpc _1 and 
redshifts z <= 1. This model will be useful for analysing 
mea urements of peculiar veloc ities at high and low redshifts 
(e.g. lHudson fc Turn bull 2012) and for studying the redshift 
space clustering signal in galaxy redshift surveys. This im- 
proved model has already been im plemented in the an alysis 
of the 6dF galaxy redshift survey i|Beutler et a/.ll2012l ). 
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